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In the Hamiltonian formulation of General Relativity the energy associated to an asymptotically 
flat space-time with metric g M „ is related to the Hamiltonian Hgr by E = Honig^v] — Hgb.\?1iiv\, 
where the subtraction of the flat-space contribution is necessary to get rid of an otherwise divergent 
boundary term. This classic result indicates that the energy associated to flat space does not 
gravitate. We apply the same principle to study the effect of zero-point fluctuations of quantum fields 
in cosmology, proposing that their contribution to the cosmic expansion is obtained computing the 
vacuum energy of quantum fields in a FRW space-time with Hubble parameter H(t) and subtracting 
from it the flat-space contribution. Then the term proportional to A* (where A c is the UV cutoff) 
cancels and the remaining (bare) value of the vacuum energy density is proportional to A 2 H 2 (t). 
After renormalization, this produces a renormalized vacuum energy density ~ M 2 H 2 (t), where 
M is the scale where quantum gravity sets is, so for M of order of the Planck mass a vacuum 
energy density of the order of the critical density can be obtained without any fine tuning. The 
counterterms can be chosen so that the renormalized energy density and pressure satisfy p = wp, 
with w a parameter that can be fixed by comparison to the observed value, so in particular one can 
chose w = — 1. An energy density evolving in time as H 2 (t) is however observationally excluded 
as an explanation for the dominant dark energy component which is responsible for the observed 
acceleration of the universe. We rather propose that zero-point vacuum fluctuations provide a new 
subdominant "dark" contribution to the cosmic expansion that, for a UV scale M slightly smaller 
than the Planck mass, is consistent with existing limits and potentially detectable. 



I. INTRODUCTION 

The cosmological constant problem is a basic issue of 
modern cosmology (see e.g. [1-4]). One aspect of the 
problem is to understand what is the effect of the zero- 
point fluctuations of quantum fields on the cosmic ex- 
pansion. Zero-point quantum fluctuations seem to give a 
contribution to the vacuum energy density of order A 4 , 
where A c is the UV cutoff. Even for a cutoff as low 
as A c ~ 1 TeV, corresponding to scales where quantum 
field theory is well tested, this exceeds by many orders 
of magnitude the value of the critical density of the uni- 
verse, which is of order (10 _3 eV) 4 . This aspect of the 
problem has by now a long history, which even goes back 
to works of Nernst and of Pauli in the 1920s [2]. 

It is clear that this A 4 contribution should somehow 
be subtracted. Renormalization in quantum field theory 
(QFT) offers a logically viable possibility. One should not 
forget, in fact, that this cutoff-dependent contribution, 
as any similar quantity in QFT, is just a bare quantity 
rather than a physical observable. Working for instance 
with a cutoff A c in momentum space one finds that, in 
Minkowski spacetime, the bare, cut-off dependent, con- 
tribution to the vacuum energy density is 



where p vac ~ (10~ 3 eV) 4 is independent of the cutoff and 
equal to the observed value of the vacuum energy density 
(assuming that vacuum energy is indeed responsible for 
the observed acceleration of the universe). The physical, 
renormalized, vacuum energy density 



PbarcW^eJ i Pcount (A c ) 



(3) 



Pbare(A c ) = cA 4 . 



(1) 



with c some constant that depends on the number and 
type of fields of the theory. The standard procedure in 
QFT is to add to it a cutoff-dependent counterterm, cho- 
sen so that it cancels the divergent part and leaves us 
with the desired finite part (for an elementary discussion 
of this point in the cosmological constant case, see e.g. 
the textbook [5]), i.e. 



Pcount (A c ) = -cA 4 + p v 



(2) 



is therefore by definition independent of the cutoff and 
equal to the observed value p va c- In principle this proce- 
dure is not different from what is done when one renor- 
malizes, say, the electron mass or the electron charge 
in QED, where again with suitable counterterms one re- 
moves the divergent, cut-off dependent, parts and fixes 
the finite parts so that they agree with the observed val- 
ues (in this sense the statement, often heard, that QFT 
gives a wrong prediction for the cosmological constant, is 
not correct. Strictly speaking QFT makes no prediction 
for the cosmological constant, just as it does not predict 
the electron mass nor the fine structure constant). 

Still, for the vacuum energy this procedure is not really 
satisfying. The source of uneasiness is partly due to the 
fact that the counterterm must be fine tuned to a huge 
precision in order to cancel the A 4 term and leave us with 
a much smaller result. In fact, if in eq. (1) we take A c 
at least larger than a few TeV, where quantum field the- 
ory has been successfully tested, Pbare is at least of order 
(TeV) 4 = (10 12 eV) 4 , so p count in eq. (3) must be fine 
tuned so that it cancels something of order (10 12 eV) 4 , 
leaving a result of order (10~ 3 eV) 4 . This fine tuning be- 
comes even much worse if one dares to take A c of the 
order of the Planck mass Mpi ~ 10 19 GeV. The fact that 
the counterterm must be tuned to such a precision creates 
a naturalness problem. Here, however, one might argue 
that neither the bare term nor the counterterm have any 



2 



physical meaning and only their sum is physical, so this 
fine tuning is different from an unplausible cancellation 
between observable quantities. The same kind of can- 
cellation appears, for instance, when one computes the 
Casimir effect. The crucial point, however, is that this 
renormalization procedure leaves us with no clue as to the 
physical value that emerges from this cancellation, so it 
gives no explanation of why the energy density associated 
to the cosmological constant appears to have just a value 
of the order of the critical density of the universe at the 
present epoch. 

The point of view that we develop in this paper is that, 
even if renormalization must be an ingredient for under- 
standing the physical effects of vacuum fluctuations, it is 
not the end of the story. The Casimir effect mentioned 
above gives indeed a first hint of what could be the miss- 
ing ingredient for a correct treatment of vacuum energies 
in cosmology. In the Casimir effect the quantity that 
gives rise to observable effects, which have indeed been 
detected experimentally, is the difference between the 
vacuum energy in a given geometry (e.g., for the electro- 
magnetic field, between two parallel conducting plates) 
and the vacuum energy computed in a reference geom- 
etry, which is just flat space-time in an infinite volume. 
Both terms are separately divergent as A^, but their dif- 
ference is finite and observable. This can suggest that, 
to obtain the physical effect of the vacuum energy den- 
sity on the expansion of the universe, one should similarly 
compute the vacuum energy density in a FRW space-time 
and subtract from it the value of a reference space-time, 
which is naturally taken as Minkowski space. 

A possible objection to this procedure could be that 
General Relativity requires that any form of energy 
should be a source for the gravitational field, which seems 
to imply that even the vacuum energy associated to flat 
space should contribute. A more careful look at the for- 
malism of GR shows however that the issue is not so 
simple and that, in a sense, the subtraction that we are 
advocating is in fact already part of the standard tenets 
of classical GR. To define carefully the energy associated 
to a field configuration in GR, it is convenient to use the 
Hamiltonian formulation, which goes back to the classic 
paper by Arnowitt, Deser and Misner [6] (ADM). As we 
will review in more detail in Section II, in order to de- 
fine the Hamiltonian of GR one must at first work in a 
finite three-dimensional volume, and then the Hamilto- 
nian takes the form 

HgR = -ffbulk + ^boundary , (4) 

where -ffbuik is given by an integral over the three- 
dimensional finite spatial volume at fixed time, and 
-^boundary by an integral over its (two-dimensional) 
boundary. At this point one would like to define the 
energy of a classical field configuration as the value of 
this Hamiltonian, evaluated on the classical solution, but 
one finds both a surprise and a difficulty. The "surprise" 
(which actually is just a consequence of the invariance 
under diffcomorphisms) is that if bulk, evaluated on any 



classical solution of the equations of motion, vanishes, so 
the whole contribution comes from the boundary term. 
The difficulty is that the boundary term, evaluated on 
any asymptotically flat metric (including flat space-time) 
diverges when the boundary is taken to infinity. The so- 
lution proposed by ADM is to subtract from this bound- 
ary term, evaluated on an asymptotically flat space-time 
with metric g^ v , the same term computed in flat space- 
time. The corresponding energy (or mass) is finite and 
is known as the ADM mass, and provides the standard 
definition of mass in GR. For instance, when applied to 
the Schwarzschild space-time, the ADM mass computed 
in this way turns out to be equal to the mass M that 
appears in the Schwarzschild metric. 

The ADM prescription can be summarized by saying 
that, in GR, the energy E associated to an asymptotically 
flat space-time with metric can be obtained from the 
Hamiltonian i?GR by 

E = H GR [g^] - H GR [rj^] , (5) 

where rj^ v is the flat metric. Even if the context in which 
this formula is valid, namely asymptotically flat space- 
times, is different from the cosmological context in which 
we are interested here, still eq. (5) suffices to make the 
point that the intuitive idea that GR requires that any 
form of energy acts as a source for the gravitational field 
is not really correct. Equation (5) tells us that the energy 
associated to Minkowski space does not gravitate. 

The idea of this paper is to generalize eq. (5) to the 
case of zero-point quantum fluctuation in curved space, 
proposing that the effect of zero-point fluctuations on 
the cosmic expansion should be obtained by computing 
the vacuum energy of quantum fields in a FRW space- 
time with Hubble parameter H (t) and subtracting from 
it the flat-space contribution. Computed in this way, the 
physical effect of zero-point fluctuations on the cosmic 
expansion can be seen as a sort of "cosmological Casimir 
effect": while in the standard Casimir effect one com- 
putes the vacuum energy in a given geometry (say, for 
the electromagnetic field, between two infinite parallel 
conducting planes) and subtracts from it the value com- 
puted in a reference geometry (flat space-time in an in- 
finite volume), here we compute the vacuum energies of 
the various fields in a given curved space-time, e.g. in 
FRW space-time, and we subtract from it the value com- 
puted in a reference space-time, i.e. Minkowski. 

A possible objection to the idea of applying eq. (5) to 
vacuum fluctuations could be that one might think that, 
after all, the structure of UV divergences is determined 
by local properties of the theory, while whether a space- 
time is Minkowski is a global question. However one 
should keep in mind that, even if T^ u {x) is a local quan- 
tity, the vacuum expectation value (0|T^„|0) in a curved 
background involves the global aspects of the space-time 
in which it is computed. This is due to the fact that 
(0|T Miy |0) requires a definition of the vacuum state |0). 
The vacuum is defined from the condition <2k|0) = 0, 
where the annihilation operators ak are defined with re- 
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spect to a set of mode functions (f>u_(t). The mode func- 
tions, in turn, are obtained solving a wave equation over 
the whole space-time, and therefore are sensitive to global 
properties of the space-time itself. In particular, in a 
FRW background the mode functions depend on the scale 
factor so their time derivatives (which enter in the com- 
putation of (0|T M „|0)) depend on the Hubble parameter 
H(t). As we will review below, this results in the fact 
that the quadratic divergence in (0|T M „|0) depends on 
the expansion rate H{t) of the FRW background. 

The paper is organized as follows. In Section II we 
recall how eq. (5) is derived in the ADM formalism. 
The reader familiar with the subject, or not interested 
in the derivation, might simply wish to move directly 
to Section III, where we apply this classical subtrac- 
tion, together with standard renormalization theory, to 
zero-point vacuum fluctuations. Some cosmological con- 
sequences of our proposal are discussed in Section IV, 
while Section V contains our conclusions. 

II. SUBTRACTIONS IN CLASSICAL GR: THE 
ADM MASS 

In this section we briefly discuss how the Arnowitt- 
Deser-Misner (ADM) mass is defined in GR [6] (we follow 
the very clear discussion of the textbook [7]). We begin 
by recalling that, in a finite four-dimensional volume V 
with boundary dV, the gravitational action is (setting 
c = 1 and using the signature r]^ = (— , +, +, +)) 

5 grav = ^J * xV =g R+j; ^ J d 3 y , 

(6) 

where g^ v is the four-dimensional metric, hij is the metric 
induced on the boundary dV, h = det hij, K is the trace 
of the extrinsic curvature of the boundary, y l are the 
coordinates of the boundary, and e = +1 on the regions 
of the boundary where dV is time- like and e = — 1 where 
dV is space-like. The first term is the usual Einstein- 
Hilbert action, while the second is a boundary term which 
is necessary to obtain a well-defined variational principle. 

To pass to the Hamiltonian formalism one performs 
the 3 + 1 decomposition of the metric, 

ds 2 = -a 2 dt 2 + hij (dx 1 + P l dt){dx ] + /3 j dt) , (7) 

where a and ft are the lapse function and shift vec- 
tor, respectively, and hij is the induced metric on the 
3-dimensional spatial hypersurfaces. One defines as usual 
the conjugate momentum tt ij = dC/dhij, where C is the 
"volume part" of the Lagrangian density, and the Hamil- 
tonian H is then the volume integral of the Hamiltonian 
density % — w^hij — C. The explicit computation (see 
e.g. Section 4.2.6 of [7]) gives 

(WnG)H = [ d 3 xVh(-aC ~-2^C') (8) 

-2 f d 2 9^ [ok - far-jiK* 3 - Kh 11 )} , 
Js t 



where S t denotes the three-dimensional spatial hyper- 
surfaces at fixed t and ft (with coordinates 6 l ) is the 
intersection of Et with dV. In an asymptotically flat 
space-time ft is just a 2-sphere at large radius r = R 
and fixed t; Co and C l depend only on hij and on its 
derivatives, but not on a and ft, while k in eq. (8) is 
the trace of the extrinsic curvature of S t , and a is the 
determinant of the two-dimensional induced metric. 

One would like to define the energy of a classical field 
configuration, i.e. of a solution of the equations of mo- 
tion, as the value of this Hamiltonian on the solution. 
Performing the variation with respect to a and ft one 
obtains the constraint equations Co = and C l = 0. 
Therefore, on a classical solution, the volume term in 
eq. (8) vanishes, and only the boundary term contributes. 
The ADM mass is defined by setting a = 1 and ft = 
after performing the variation (corresponding to the fact 
that energy is associated to asymptotic time transla- 
tions; setting a = and ft = 1 one rather obtains 
the ADM momentum F| DM ), so only the term propor- 
tional to k contributes to the mass. However, even for 
Minkowski space, this boundary term diverges. In fact, 
for an asymptotically flat space-time we can take dV to 
be a three-dimensional cylinder made of the two three- 
dimensional time-like hypersurfaces {t = t\ , r < R} and 
{t = t 2 ,r < R} (the "faces" of the three-dimensional 
cylinder) and of the space-like hypersurface {r = R, t\ < 
t < <2}- Let us denote by Kq the extrinsic curvature of 
dV computed with a flat Minkowski metric. The faces 
at t = t\ and ti have Kq — 0; however, on the surface 
{r = R, t\ < t < £2} the extrinsic curvature is Kq = 2/R, 
while = R 2 sm 2 9 and e = +1, so the boundary 

term in eq. (7) is [7] 

/ d 3 y ^/\h\eK = 87r(i 2 - tx)R , (9) 

JdV 

and diverges both when R — > 00 and when (i2 — ii) 00. 
As a consequence, also the boundary term proportional 
to k in eq. (8) diverges, already for the Minkowski met- 
ric, and then of course also for generic asymptotically flat 
space-times. The ADM prescription is then to replace K 
in eq. (6) by (K — K ), i.e. to subtract from the trace of 
the extrinsic curvature computed with the desired met- 
ric, the value computed in flat space. Correspondingly, 
the trace k of the two-dimensional extrinsic curvature 
in eq. (8) becomes (k — k ), and the ADM mass of an 
asymptotically flat space-time is defined as [6] 

M ADM = --^— lim / d 2 9^(k-k ). (10) 

07TG St^oo J g 

If for instance one applies this definition to the 
Schwarzschild space-time, one finds that Madm is equal 
to the mass M which appears in the Schwarzschild met- 
ric. 

What we learn from this is that, in classical GR, the 
mass or the energy that acts as the source of curvature 
of space-time, such as for instance the mass M that en- 
ters in the Schwarzschild metric, can be obtained from a 
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Hamiltonian treatment only after subtracting a flat-space 
contribution that need not be zero, and is in fact even 
divergent, but still does not act as a source of curvature. 



III. APPLICATION TO ZERO-POINT 
ENERGIES 

It is natural to apply the same principle to zero-point 
quantum fluctuations. In particular, in a cosmological 
setting, we propose that the zero-point energy density 
and pressure that contribute to the cosmological expan- 
sion are obtained by computing the energy density and 
pressure due to zero-point quantum fluctuations in a 
FRW metric with Hubble parameter H(t) = a/a, and 
subtracting from it the flat-space contribution. 

We consider first the contribution of a real massless 
scalar field. In a FRW background the mode expansion 
of the field is 



(*) = / — 
J (2ir 



d 3 k 



(27r) 3 V2fc 



a k ^(t)e ik ' x + 4^(t)e" 



tkx 



(11) 

where k is the comoving momentum, and <j>k{t) satis- 
fies the massless Klein- Gordon equation in a FRW back- 
ground, 



b'l + 2-<j>' k + k 2 <f> k = Q, 



(12) 



where the prime is the derivative with respect to confor- 
mal time rj. Writing cf) k — ip k /a this equation is reduced 
to the form 

^'+(fc 2 -^)^ = 0. (13) 

In a De Sitter background we have a(t]) = —\/{Hr]) so 
a" j a = 2/rj 2 , while in a matter-dominated (MD) epoch 
a ~ rj 2 and therefore again a" /a = 2/rj 2 . Thus, in both 
cases the positive frequency solution of eq. (13) is 



Mv) = i - 



kr] 



-ikrj 



and therefore ^>k(v) 1S given by 
4>k(ri) 



1 \ — ) e~ ikl ' 



a(rj) \ kr\ 



(14) 



(15) 



with a(rj) the scale factor of De Sitter or MD epoch, 
respectively. In contrast, during a radiation-dominated 
(RD) epoch a ~ rj, so a" /a = 0. Then ip k (r]) = cT^ 
and 



(16) 



Using the energy-momentum tensor of a minimally cou- 
pled massless scalar field, 



(17) 



and setting g^ v = (— 1, a 2 Sij), a simple computation [8, 
9], reviewed in Appendix A, shows that the off-diagonal 
elements of (0|T M „|0) vanish, while the vacuum energy 
density and pressure are given by 



d 3 k 



(2ir) 3 2k 



\M 2 + %\M 2 

a z 



3a 2 



(18) 
(19) 



where the dot denotes the derivative with respect to cos- 
mic time t. As they stand, these expressions must still be 
regularized, and we regularize them by putting a cutoff in 
momentum space. Recall that the comoving momentum 
k is just a label of the Fourier mode under consideration, 
while the physical momentum of the mode is given by 
k/a. We expect that quantum gravity enters the game 
when the physical momenta exceed the Planck scale, and 
we therefore put a time-independent cutoff A c over phys- 
ical, rather than comoving, momenta. In terms of co- 
moving momentum k this means k < a(t)A c . Since the 
modes (f)k depend only on k — |k|, the angular integrals 
are trivially performed, and we finally get 



aA c 



p=-y o d kk^ k \ 2 + k ^ 2 



P = 



8tt 2 



dk k 



(20) 



(21) 



A. Vacuum fluctuations in De Sitter space 

1. Vacuum energy density 

We compute first these expressions in De Sitter space. 
We therefore plug eq. (15), with 0(77) = —l/{Hrf), into 
eqs. (20) and (21). For the energy density the result is 



p{K) 



At 



+ 



dk k 
H 2 A 2 



16tt 2 16tt 2 



(22) 



The first term is the well known flat-space result, propor- 
tional to the fourth power of the cutoff, while the term 
quadratic in A c is the correction due to the expansion of 
the universe. The appearance of a term <~ H 2 A 2 C , in a 
theory with two scales, the UV cutoff A c and the Hub- 
ble scale H, can also be understood using rather general 
arguments [10]. 

According to our proposal we now subtract the flat- 
space contribution, which is simply the term ~ A* in 
eq. (22), and we find that in De Sitter space the zero- 
point quantum fluctuations of a real scalar field have a 
(bare) energy density 



Pbare(A c ) 



g 2 A 2 
16tt 2 



(23) 
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The subscript "bare" stresses that this is still a bare, cut- 
off dependent quantity. We have eliminated the quartic 
divergence thanks to the "classical" prescription (5), but 
the result is still a bare energy density, diverging quadrat- 
ically with the cutoff. In this sense, the situation is dif- 
ferent from the usual Casimir effect, where the subtrac- 
tion of the flat-space contribution suffices to make the 
result finite. However, we can now use standard renor- 
malization theory, so the renormalized energy density is 
obtained by adding a counterterm whose divergent part is 
chosen so to cancel this quadratic divergence, and whose 
finite part p vac is in principle fixed by the observation, so 



Pcount(A c ) 



H 2 A 2 
16tt 2 



+ Pv 



(24) 



As usual in the renormalization procedure, the finite part 
/3 vac cannot be predicted. It must be fixed to the observed 
value. What we have gained, with respect to eq. (2), is 
that now we have a different understanding of what is 
a "natural" value of this finite part. If quantum gravity 
sets in at a mass scale M (so that M could be typically 
given by the Planck mass Mpi, or by the string mass), a 
natural, non fine-tuned value of p vac is given by 



const, x a- 



H 2 M 2 
16tt 2 



(25) 



where "const." is a constant 0(1), which cannot be fixed 
by naturalness arguments only, and must be determined 
by comparison with the experiment. [41] 

Observe that naturalness arguments cannot fix the sign 
of the finite part either, and we have added a factor a = 
±1 to take this fact explicitly into account. There is 
in fact no a priori reason why the renormalized vacuum 
energy should necessarily be positive. For instance, the 
vacuum energies obtained for fields in a finite volume 
from the Casimir effect can be either positive or negative, 
depending e.g. on the type of field and on the geometry 
considered. 



2. Equation of state and general covariance 



Repeating the same computation for the pressure we 



get 



p(A c 



^ 2 l 



dk k 



4tt 2 J \ 3a 4 6a 2 

A 4 H 2 A 2 



48tt 2 48tt 2 



(26) 



The term ~ A 4 in cq. (26) was already computed in [11] 
and is just the result in Minkowski space, that we sub- 
tract, so we end up with 



H 2 A 2 
' 48tt 2 



Pbare(A c ) 



(27) 



Observe that Pbarc(A c ) = -(l/3)/9 b arc(A c ). However, it 
would be incorrect to conclude that vacuum fluctuations 



in De Sitter space satisfy the equation of state p = wp 
with w = —1/3. The point is that this relation only 
holds for the bare quantities, and not necessarily for the 
renormalized ones. The physical, renormalized, pressure 
is obtained by adding a counterterm 



Pcount(A c ) = + 



2A2 



H 2 A^ 
48tt 2 



(28) 



Observe that regularizing the theory with a cutoff over 
spatial momenta, as we have done, breaks explicitly 
Lorentz invariance in Minkowski space, since the no- 
tion of maximum value of spatial momenta is not in- 
variant under boosts. In a generic FRW background, of 
course, Lorentz transformations are not a symmetry of 
the theory, since the metric depends explicitly on time, 
and the guiding principle is rather general covariance, 
which again is broken by a cutoff over spatial momenta. 
However, even if the regularization breaks general covari- 
ance, in De Sitter space a generally covariant result can 
still be obtained in the end for the physical, renormal- 
ized, vacuum expectation value of the energy-momentum 
tensor, just by choosing the finite parts of the countert- 
erms such that p vac = —p vac - Then the vacuum expecta- 
tion value of the renormalized energy-momentum tensor 
= dia,g(—p,p,p,p) becomes 



(0\T»\0) = const, x cr 



H 2 M 2 
16tt 2 



(29) 



or, lowering the upper index with the metric .g M „, 



(0|T^|0) = const, x cr 



H 2 M 2 
16tt 2 



(30) 



Since in De Sitter space H is constant we see that, with 
the choice p vac — —p vac , (0\T^\0) is just given by a nu- 
merical constant times g^ v , and it is therefore covari- 
antly conserved. In the language of the effective action 
for gravity, which is obtained by treating the metric g^ u 
as a classical background and integrating over the mat- 
ter degrees of freedom (see e.g. refs. [12-14]), the vac- 
uum expectation value of the energy-momentum tensor 
is given by the functional derivative (2/y^g)S/Sg^ u of 



the effective action. Then a contribution such as that 
given in eq. (30) can be obtained by taking the functional 
derivative of the volume term in the effective action (see 
rcf. [15] for a discussion of the equation of state that can 
be obtained from the various contributions to the effec- 
tive action). 

The choice p vac = — Pvac will therefore be assumed in 
the following, for De Sitter space-time. Observe that a 
covariant result for the renormalized value of (0|T M „|0) is 
obtained with a counterterm that is not covariant, i.e. is 
not proportional to g^, since p CO unt(A c ) ^ -p CO unt(A c ), 
but again this is a consequence of the fact that our reg- 
ularization is not covariant. 
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3. Contribution of higher-spin fields and massive particles 

A similar conclusion about the "natural" value of the 
energy density of vacuum fluctuations holds if we add the 
contribution of fields with different spin. In fact, massless 
fermions and gauge bosons have a conformally invariant 
action and then, since the FRW metric is conformally 
equivalent to flat Minkoswki space, their vacuum energy 
in a FRW space-time is the same as in Minkowski space. 
Therefore, with our subtraction (5), they do not con- 
tribute to pvac and pvac- This is completely equivalent 
to the well known fact that vacuum fluctuations of mass- 
less fermions and gauge bosons are not amplified during 
inflation. 

The contribution of massive fermions to vacuum fluc- 
tuations in De Sitter space, as well as the generalization 
of eq. (22) to massive bosons, can be readily computed. 
For massive bosons, eq. (22) becomes [9, 16] 

/>£lA c 

1 



1 



Pb(A c ) = ^ 



dk k 



kuj k 

~3~ 



2a 3 Wfc 



where we have neglected terms of order H 4 . 
convergent in the UV; lju is defined as 



\J m 2 



(k/a) 



(31) 
and terms 



(32) 



so in the massless limit ook — > k/a and we recover eq. (22) 
[42]. In this massive case our prescription amounts to 
subtracting the flat-space contribution given by the term 
fcwfe/a 3 in brackets. The remaining, quadratically diver- 
gent term, leads again to eq. (25), times a correction 
[1 + C(m 2 /M 2 )], and to a logarithmically divergent term, 
which after renormalization produces a contribution to 
Pvac proportional to H 2 m 2 log(M/m), and therefore sub- 
leading for m <C M. 

The contribution of a massive Majorana spinor field 
to the vacuum energy density in De Sitter space is in- 
stead [16] 



Pf(A c ) = 



2tt 2 



aA c 



dk k 



kuj k 



+ 



m 2 H 2 k 



8a 3 u 3 



(33) 



The first term in brackets gives the usual negative con- 
tribution to vacuum energy in flat space due to fermions. 
Observe that in a supersymmetric model, where to each 
Majorana spinor is associated a complex scalar field, 
and therefore two real scalar fields, the contribution 
(—kcok/a 3 ) in eq. (33) cancels exactly the bosonic con- 
tribution (+kuj k /a 3 ) in eq. (31), giving the usual can- 
cellation of vacuum energy for a supersymmetric theory 
in Minkowski space. In a realistic model with super- 
symmetry broken at a scale A susy , this cancellation is 
however only partial and leaves the usual result of order 
Af usy . In our approach, instead, the quartic divergence 
are eliminated exactly by the prescription of subtracting 
the flat-space contribution. 

The term that remains in eq. (33), after subtraction 
of the flat-space contribution, gives a divergent contri- 
bution equal to (m 2 H 2 /16tt 2 ) log(A c /m) (consistently 



with the fact that, for m = 0, it must vanish because 
of conformal invariance). After renormalization, this 
gives a contribution to the vacuum energy which is of 
order (m 2 H 2 / 16tt 2 ) \og(M / m) , and which therefore for 
m <C M, is negligible with respect to the bosonic contri- 
bution (25). It is also interesting to observe that, even in 
a theory with exact supersymmetry, the cancellation be- 
tween fermionic and bosonic divergences only takes place 
at the level of the quartic divergence. There is no contri- 
bution proportional to A 2 from the fermionic sector, and 
therefore the whole contribution proportional to H 2 A 2 
comes from the bosonic sector. This means that eq. (25) 
holds even in a theory with exact or broken supersym- 
metry. 

In contrast, gravitons give a contribution to the bare 
energy density proportional to A 2 H 2 . In fact, in a FRW 
space-time each of the two helicity modes h a ,k, with 
a = {+, x} satisfies separately the same wave equation 
as eq. (12) with replaced by h a> k. 



K,k + 2-h' aik + k 2 h a , k 



0. 



(34) 



So each of the two helicity modes gives the same contri- 
bution to p vac and p vac as a minimally coupled massless 
scalar field. Therefore, in a theory with n s minimally 
coupled elementary scalar fields plus the two degrees of 
freedom for the graviton, the natural value for the energy 
density, eq. (25), is of order (n s + 2)H 2 M 2 /(16tt 2 ). In 
particular, even in a theory with no fundamental scalar 
field, a contribution of order H 2 M 2 /(8tt 2 ) comes anyhow 
from gravitons. 



B. Vacuum fluctuations during RD and MD 

It is straightforward to repeat the same calculation for 
a radiation-dominated (RD) and for a matter-dominated 
(MD) era. For the RD epoch we use the modes (16). 
Then, recalling that dt = adrj, 



4>k = -<t>' k = — ~o 



a 

ik-\ 

a 



ikr} 



(35) 



where as usual the dot denotes the derivative with respect 
to cosmic time t and the prime the derivative with respect 
to r\. Using a 1 / a = a = aH(t), we get 



P(A C ) 



1 

4^2 



dk k 



H 2 (t) 
2a 2 



K , H 2 (t)A 2 c 



16tt 2 



16tt 2 



(36) 



so the energy density turns out to be identical to eq. (22), 
except that now H is replaced by H (t). For the pressure 
we get 



P(A C 



K , H 2 (t)A 2 c 



48tt 2 16tt 2 ' 



(37) 
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Thus, once removed the Minkowski term, we remain with 



Pbare(A c ) 



16tt 2 



(38) 



while for the pressure we get pbare(A c ) = /9bare(A c ). Sim- 
ilarly to what we have discussed in the dc Sitter case, this 
means that the natural value of the renormalized energy 
density is 



Pvac(i) = const, x a 



H 2 {t)M 2 
16tt 2 



(39) 



where M is the scale where quantum gravity sets in, 
"const." is a numerical constant 0(1), and a = ±1. On 
the other hand, the relation Pbare(A c ) = Pbare(A c ) does 
not imply that the renormalized energy density and the 
renormalized pressure satisfy an equation of state with 
w = +1. As in the de Sitter case, the finite part in 
the counterterm for the pressure can be chosen so to re- 
produce any observed value of w, in particular the value 
w = —1. The issue of general covariance is however now 
more subtle, since even the choice w — — 1 now leads to 



(0\T^\0) = const, x cr 



H 2 (t)M 2 
16tt 2 



(-9^) (40) 



which, because of the time dependence of H(t), is no 
longer covariantly conserved. However, general covari- 
ance actually only requires that the total energy momen- 
tum tensor, including that of matter, radiation, etc., be 
covariantly conserved. The fact that energy-momentum 
tensor associated to vacuum fluctuations is not separately 
conserved means that there must be an energy exchange 
between vacuum fluctuations and other energy sources 
such as radiation or matter, so that we are actually deal- 
ing with an interacting dark energy model. We will come 
back to this issue in Section IV. 

For MD we use the modes given in eq. (15), with a(rj) <~ 
rj 2 . For the energy density we find 



P(A C ) = 



1 

4^ 



aA c 



dk k 



H 2 {t) 9H A (t) 



2a 2 



32k 2 



(41) 

So, both for MD and for RD, the first two terms are the 
same as in De Sitter, eq. (22), except that H becomes 
H(t). Again, the first term is the Minkowski contribu- 
tion, that we subtract. The term ~ H A diverges only 
logarithmically with the UV cutoff (and also require an 
IR cutoff, which for a light scalar field is provided by its 
mass, while for a strictly massless field could be taken 
equal to H). Thus, for the physical, renormalized, en- 
ergy density we find again that the natural value is given 
by eq. (39), since the term <~ H A in eq. (41) produces a 
subleading term of order H A (t) \og(M / H(t)). 

Recalling that (in units K = c = 1) Newton's constant 
is G = 1/Mpj, and that the critical density at time t is 



Pc(t) 



3H 2 (t) 
8irG 



8tt 



H 2 {t)M 2 x , 



(42) 



we can express the above result by saying that the "nat- 
ural" value suggested by QFT is 



Pvac(i) = const, x a 



1 



M 



Pc(t) , (43) 



with const. = 0(1) and a = ±1. We therefore find that, 
both during RD and MD, the renormalized energy den- 
sity due to zero-point quantum fluctuations is not a con- 
stant, and therefore does not contribute to the cosmolog- 
ical constant. Rather, it is a fixed fraction of the critical 
density p c {t). 

As before, each helicity mode of a graviton contributes 
as a minimally coupled scalar field, while gauge bosons 
do not contribute and the contribution of light fermions 
is suppressed by a factor {m / M) 2 \og(M / m) . Thus, in 
a theory with n s fundamental minimally coupled scalar 
fields plus the two helicity modes of the graviton we have, 
for either De Sitter, RD or MD, 



p z (t) = £lzPc(t) ■ 



(44) 



where the subscript Z stands for "zero-point quantum 
fluctuations" , and the natural value of Viz is 



ft; 



, K + 2) / m y 



(45) 



where M is the scale where quantum gravity sets in, e.g. 
Mpi itself or the string scale. As we have seen in eqs. (31) 
and (33), bosons or fermions with mass m not far from 
the quantum gravity scale M would give corrections of or- 
der [1 + C(m 2 /M 2 )] to this estimate (which could become 
numerically important if there were, e.g. many fermions 
at a scale, such as the GUT scale, not far from the Planck 
mass). The computation of the pressure during MD gives 



P(A C ) = 



1 

I^ 2 



f k 2 | H 2 (t) i 9H 4 (t) 



dkk {3a~ A+ 3a 2 



+ 



32fc 2 



(46) 

Again, the first term is the Minkowski contribution, that 
we subtract, and we neglect the term ~ H A . Thus, we 
finally find 



Pbare(Ac) = 



H\t)K 2 c 
24tt 2 



(47) 



and therefore, during MD, Pbaro = +(2/3)/5barej but again 
the renormalized energy and pressure satisfy p — wp with 
w determined by the observation. [43] 

A technical point that deserves some comment is the 
choice of the modes given in eqs. (15) and (16). These 
modes are particularly natural since in the UV limit they 
reduce to positive-frequency plane waves in flat space. 
However, the choice of the modes is equivalent to the 
choice of a particular vacuum state, and the most general 
possibility is a superposition of positive- and negative- 
frequency modes with Bogoliubov coefficients ak and (3k ■ 
As we show in appendix B, for a generic choice of vacuum 
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the dependence of the natural value of p vac on H 2 (t)M 2 
is not altered, while the numerical coefficient in front of 
it can change. 

A conceptually interesting aspect of the result (39) is 
that it appears to involve a mixing of ultraviolet and in- 
frared physics, since it depends both on the UV scale M, 
and on the horizon size if -1 , which represents the "size 
of the box" , and therefore plays the role of an IR cutoff. 
This is an interesting result by itself, since in quantum 
field theory we are rather used to the fact that widely sep- 
arated energy scales decouple. Observes that this UV-IR 
mixing comes out only because our classical subtraction 
procedure based on eq. (5) eliminates the troublesome 
term diverging as A^, which is instead a purely UV term. 
As we already discussed in the Introduction, the origin of 
this UV-IR mixing can be traced to the fact that, even if 
(x) is a local quantity, the vacuum expectation value 
(0|T M1/ |0) is sensitive to H(t) through the definition of the 
vacuum state |0). In fact, the vacuum is defined from the 
condition ak|0) = 0, where the annihilation operators ak 
arc defined with respect to a set of mode functions 0k (i); 
these mode functions are obtained solving a wave equa- 
tion over the whole space-time, and therefore are sensi- 
tive to the time evolution of the scale factor and, more 
generally, to the overall geometry of space-time. 



IV. COSMOLOGICAL IMPLICATIONS 

A. Vacuum fluctuations and the dominant 
component of dark energy 

The results of the previous section show that, with the 
subtraction that we advocate, zero-point vacuum fluc- 
tuations do not contribute to the cosmological constant 
since their energy density, given in eqs. (44) and (45), is 
not a constant, but rather a fixed fraction of the critical 
energy at any epoch. 

The first question to be addressed is whether a vac- 
uum energy density with such a time behavior could 
be identified with the dark energy component Qa which 
is responsible for the observed acceleration of the uni- 
verse. If this were the case, since we have found that 
the vacuum energy density scales as H 2 (t), we would 
have flA(t) = &,AH 2 (t)/H 2 where we follow the stan- 
dard use of denoting by J7a the value of fl\(t) at the 
present time t = io- Such a model has been compared to 
CMB+BAO+SN data in rcf. [17], where it is found that 
it is ruled out at a high significance level. 

Formally the result of ref. [17] comes from the fact that 
the x 2 obtained fitting this model to the usual estimators 
for CMB, for BAO and to SNIa turns out to be unaccept- 
ably high. However, independently of technical details, 
it is easy to understand physically why such a model for 
dark energy is not viable. Consider in fact a spatially 
flat model with vacuum energy density PA(t) evolving as 
H 2 (t), with equation of state w\ = — 1, and matter den- 
sity PAi(t), with Q,m + &a = 1, in the recent universe 



where we can neglect fi#. The total energy-momentum 
conservation is 



PM + PA + 3Hp M = , 



(48) 



and cannot be split into two separate conservation equa- 
tions for pa and pm, Pm + ^>HpM = and pa = 0, 
since the second equation is obviously incompatible with 
PA(t) ~ H 2 (t), in the recent universe where the Hubble 
parameter is certainly not a constant. This means that 
energy must be transferred between pA(t) and pM(t) in 
order to obtain the behavior pA(t) ~ H 2 (t), so we have 
an interacting dark energy model. The Friedmann equa- 
tion in this model reads 



H 2 (t) 
Hi 



and therefore 



and 



H 2 {t) _ 1 



h 2 i - n A 



n M (i) 



n A {t) = 



_TLa_ 

1-flA 



(49) 



(50) 



(51) 



In this model, therefore, the time evolution of dark en- 
ergy density is the same as that of matter. This is 
clearly incompatible with the existing cosmological ob- 
servations, that rather indicate that in the recent epoch 
&A(t) remained constant, at least to a first approxima- 
tion, while the matter density in the concordance ACDM 
model evolves in a way that cannot differ too much 
from Om(£) ~ l/a 3 (t). More quantitatively, combin- 
ing the Friedmann equation (50) with the total energy- 
momentum conservation equation (48) it can be easily 
found [17] that in this model the dependence of 17m (i) 
on the scale factor a(t) is 



n M (t) 



n 



M 



a 



3-3Q A 



(52) 



where we normalize the scale factor a(t) so that a (to) = 1 
at the present epoch to (we will see explicitly in Sect. IV C 
how to derive this result in a similar setting). Phys- 
ically this result expresses the fact that in this model 
Pa(£) ~ H 2 (t) decreases with time, rather than remain- 
ing constant as it would do in isolation, and therefore 
part of its energy density must be transferred to matter, 
through the energy-conservation equation (48). So, mat- 
ter energy density decreases more slowly that 1/a 3 . For 
flA{t), recalling that we are considering a flat model with 
ft A + Hm = 1) eq. (51) gives 



n A (t) 



0. 



t3 — 3S"2a 



(53) 
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To compare these results to observations, we recall that 
the possibility of a time evolution of the dark energy den- 
sity is usually studied in the literature by parametrizing 
it as 



n A (t) 



a 3+3w A (t) 



(54) 



It is important to stress that eq. (54) is simply an effective 
way of parametrizing the time dependence of pde in order 
to compare it with the observations, and the parameter 
i«A that appears there is equal to the w parameter that 
enters in the equation of state of dark energy only if dark 
energy is not interacting which, as we have seen, is not 
the case for a model where p\(t) <~ H 2 (t). 

Comparing eqs. (53) and (54) we see that a model with 
Q,A (i) = Q\H 2 (t) I Hq predicts that the effective parame- 
ter w\ in eq. (54) is constant in time, and equal to — Ct\. 
The limit on constant wa obtained from WMAP 7yr 
data+BAO+SN is w A = -0.980 ± 0.053 at 68% c.l. [18] 
(actually, this value does not include systematic errors in 
supernova data; including the systematics, the error on 
wa rather becomes about 0.08, see [19]), and therefore re- 
producing this value would require fi A = +0.980 ±0.053, 
and Qm = 1 — consistent with zero to a few per- 
cent! It is no wonder that such a model does not fit 
(any) other cosmological observation. Alternatively, set- 
ting Qm — 0.26 and Qa — 0.74, we get a prediction 
wa — —0.74, which is excluded at a high confidence level. 

Therefore, it appears that the interpretation of zero- 
point quantum fluctuations as the dominant dark energy 
component, responsible for the acceleration of the uni- 
verse, is not viable. In the next section we will explore 
the possibility that they could still provide a new, sub- 
dominant, "dark" component, that we will denote by pz 
to distinguish it from the dominant component pa which 
instead, according to the limits on wa, is at least approx- 
imately constant in time. We will see that, for plausible 
values of the mass scale M where quantum gravity sets 
in, the energy density pz can have a value consistent with 
existing observations, but still potentially detectable. 



B. Theoretical expectations for Qz 

The effect of pz on the cosmological expansion depends 
on the mass scale M, whose exact value can only be de- 
termined once one has a fundamental theory of quantum 
gravity. If we set M = Mp\ and we consider the Stan- 
dard Model with one Higgs field, so n s = 1, eq. (45) gives 
\Qz\ — l/(27r) ~ 0.16. However, precise numerical fac- 
tors are beyond such order-of-magnitude estimates and, 
by lowering the UV scale M, it is easy to reduce this 
number to smaller but still potentially observable values. 
For instance, in heterotic string theory the scale is rather 
given by the heterotic string mass scale Mh = gMp\, 
where g ~ 1/5 is the value of the gauge couplings at the 
string scale [20]. This would rather lead to the estimate 
\Q Z \ ~ (n s + 2)g 2 /(6n) ~2x 10 _3 (n s + 2). Lower val- 
ues of the cutoff, possibly down to the TeV scale, can 



be obtained in theories with large extra dimensions [21]. 
It is also important to observe that, as explained in the 
discussion below eq. (33), the estimate given in eq. (45) 
holds even in a theory with exact or broken supersymme- 
try, since the contribution proportional to H 2 M 2 comes 
anyhow only from the bosonic sector. 



C. Cosmological evolution equations 

We consider a fiat ACDM cosmology, with a vacuum 
energy Qa having an equation of state pa = waPa and 
we further add the energy density p z {t) given in eq. (44), 
with pz = wzpz- We have nothing to add to the problem 
of the physical origin of p A , except that in our model it is 
not due to zero-point quantum fluctuations: it has noth- 
ing to do with the quartic divergence in the vacuum en- 
ergy (which is eliminated by our ADM-like subtraction), 
nor with the quadratic divergence, which is instead the 
origin of pz- Then in this model (which could be conve- 
niently called AZCDM) 



p = Pr + Pm + PA + Pz , 
1 

P = ^Pr + WA ? A + Wz P z ■ 



(55) 
(56) 



The values wz = wa = — 1 will be assumed in the fol- 
lowing (in appendix C we discuss the case wz generic). 
From eq. (44), 



p z (t) = MzPc(t) 



^zPo \ Tf2 
H 2 



H 2 (t), (57) 



where po = 3Hfi / (8ttG) is the present value of the critical 
density, so the energy density in the dark sector is 



Pde 



Pa + Pz 



= PA 



3Q Z 
8irG 



H 2 {t) 



(58) 



Quite interestingly, this is the same form of the dark 
energy density found in refs. [22-25] from an appar- 
ently rather different approach, namely from the sugges- 
tion that the cosmological constant could evolve under 
renormalization group, after identifying their parame- 
ter v with our Qz, compare with eq. (13) of rcf. [23]. 
Observe also that their value for the parameter v is 
v = ±(l/127r)M 2 /Mp[ where M is the mass scale where 
new physics comes in. This has the same parametric de- 
pendence on (M/Mpi) as our result for Qz, and is even 
quite close numerically [44]. 

Using eq. (58), the Friedmann equation becomes 



H 2 (t) = 



8ttG 



Pr + Pm + PA + 



3Qz 
8ttG 



H 2 (t) 



, (59) 



which can be rewritten as 
Hi 



H 2 (t) = —±-[n R (t) + n M (t) + nA(t)}, (60) 
i — \i z 
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where fl^t) = pi{t)/p , (i = R,M,A). Thus, the ef- 
fect of £lz on the Friedmann equation is equivalent to 
a rescaling of the present value of the Hubble constant, 
Hq — »■ Hq/(1 — Qz) 1 ^ 2 or, cquivalently, to a rescaling of 
fli(t) (with i = R,M,A) into - fl z ). Even if 

we set wa = — 1 we have for the moment written Q\(t), 
rather than setting it to a constant, in order to allow for 
the possibility of an energy exchange between pz(t) and 
p\(t), see below. 

Zero-point fluctuations also contribute to the equation 
for the acceleration 



a 
a 



(fi + 3p) , 



(61) 



in a way which depends on the sign of pz ■ With wz = — 1 
(or more generally whenever wz < —1/3), zero-point 
fluctuations contribute to accelerating the universe if 
pz > 0, while for pz < they give a contribution that 
decelerates the expansion, and which therefore opposes 
the accelerating effect of pA- 

We next consider the energy conservation equation. 
As we already mentioned in Section IIIB, the fact that 
zero-point quantum fluctuations have an energy density 
Pz(t) ~ H 2 {t) has non-trivial consequences on the con- 
servation of energy. Consider in fact the energy conser- 
vation in a FRW background, 



p = -3H{p + p). 



(62) 



When there is no exchange of energy among the differ- 
ent components, energy conservation is satisfied sepa- 
rately for each component, so p, = —3H(pi + pi) with 
i = R, M, A, Z. If this were the case, using pz = wzPz, 
the equation 



pz = -3H(p z + pz) = -3(1 + w z )Hp z 



(63) 



would give p z (t) <~ a~ 3 ( 1+Wz \ and in particular p z (t) 
constant if w z = —1. However, we have found that 
Pz(t) ~ H 2 (t). So, energy must be transferred between 
zero-point fluctuations and other components [45]. 

In principle, there are various mechanisms by which 
vacuum fluctuations can exchange energy with ordinary 
matter. A typical example is the amplification of vacuum 
fluctuations [27, 28], or the change in a large-scale scalar 
field due to the continuous flow across the horizon of 
small-scale quantum fluctuations of the same scalar field, 
which is also at the basis of stochastic inflation [29] . If we 
assume that pz exchanges energy with pm but not with 
Pa the relevant conservation equation is (setting hereafter 
w z = -1) 



Pm + Pz = -3Hp M ■ 
Using eq. (57), we then obtain 

'n z p a \ dH 2 (t) 



Pm 



-3Hp 



M 



Hi 



dt 



(64) 



(65) 



This equation was already discussed in refs. [17, 24] in 
the context of their time-varying cosmological constant 



model and, following these papers, to solve it we compute 
the term dH 2 (t)/dt on the right-hand side by using the 
Friedmann equation (60) and we obtain (neglecting for 
simplicity Q.n(t) in the low redshift epoch in which we 
are interested here) 



Pm 



-3Hp 



M 



1-fiz 



PM ■ 



(66) 



where we used the fact that, since we are assuming that 
pz only interacts with pm, and we are furthermore as- 
suming w\ = — 1, the energy density pa evolves in isola- 
tion and satisfies pa — 0. Equation (66) can be rewritten 
as 

p M = -3(1 - n z )H PM , (67) 

and has the solution 

P mW-Pm(0)(1 + z) 3(1 - S!z) . (68) 

Of course, since we have taken wz = —1, if vacuum en- 
ergy density were non-interacting it would remain con- 
stant in time. In the case a = +1 (i.e. when Qz > 0) we 
have rather found that pz is proportional to +H 2 (t) , and 
therefore it decreases with time instead of staying con- 
stant. This means that dark energy is interacting, and 
that energy is transferred from vacuum fluctuations to 
matter, for instance with a mechanism analogous to the 
amplification of vacuum fluctuations, and as a result the 
energy density of matter must decrease slower than 1/a 3 . 
This is reflected in eq. (68), since for Viz > we find 
that pm indeed decreases slower than 1/a 3 . If Qz < 
the situation is reversed. A behavior p z (t) ~ —H 2 {t) 
means that p z (i) becomes less and less negative as time 
increases, so energy is transfered from matter to vacuum 
fluctuations, and SIm{z) decreases faster than 1/a 3 . 

In the absence of an understanding of the dynamical 
origin of the dominant dark energy term pa, it is inter- 
esting to consider also the possibility that energy could 
be exchanged also between pz and pa- We assume at 
first, for simplicity, that energy is exchanged only with 
Pa, and not with pu- Then the corresponding conserva- 
tion equation (taking wa = — 1 for simplicity) is 



PA + Pz = , 
which trivially integrates to 

P\(t) + pz(t) = constant = pa{U)) + Pzih) 
and therefore, using pz{t) = flzPoH 2 (t)/ , 

■H 2 (z) 



Qa(z) = Qa-&z 



Hi 



(69) 



(70) 



(71) 



where, as usual, on the right-hand side = Qa{z = 0). 
However, in terms of the total dark energy density de- 
fined in eq. (58), we see that in this case we simply have 
a total dark energy density pde that satisfies />de = 0, 
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while the matter energy density satisfies its usual conser- 
vation equation pm + 3HpM = 0, and also eqs. (59) and 
(61), when rewritten in terms of pde> take the standard 
ACDM form. Thus, in the end a model where pz only ex- 
changes energy with pa (and in which wz = wa — — 1) is 
indistinguishable from standard ACDM cosmology with 
just a cosmological constant. In the most general case 
in which wz + 1 —1, however, even setting w\ = — 1 
the pressure pde = Pz + PA = u'zPz + waPa is no 
longer equal to — Pde, so this model has observable de- 
viations for ACDM, and is in fact of the type called 
AXCDM [30-33]. In the following we will however re- 
strict to wz = wa = — 1. 

A more general phenomenological analysis, in which 
one takes into account the possibility that pz interacts 
both with pa and with pm, can be performed by splitting 
the conservation equation 



Pm + PA + Pz 



-3Hp 



M 



into the two equations 

PA = 
Pm = 



-(1 - a)p z , 
-'SHp M - ap z ■ 



(72) 



(73) 
(74) 



where < a < 1. Equation (64) corresponds to the 
limiting case a = 1 while eq. (69) corresponds to the 
limiting case a = 0. Rewriting these equations in terms 
of the total dark energy density Pde = PA + Pz, we get 



p DE = ap z , 
Pm = -3Hp M - ap z ; 



(75) 
(76) 



which shows that the observable consequences depend 
only on the combination apz- In other words, as long 
as wz = wa = — 1, there is no point in postulating an 
energy exchange between pz and pa, since only the frac- 
tion of pz which is exchanged with pM has observable 
consequences. In the following we will limit ourselves to 
wz = wa = — 1, and we then set a = 1. 



D. Limits on Qz from cosmological observations 

In this section we perform a first analysis of the lim- 
its that some cosmological observations impose on VLz- 
A more detailed comparison with the data will be pre- 
sented elsewhere. The fact that the energy budget of 
the universe at the present epoch is known to a precision 
of about 1% by itself does not yet constraint £lz, since a 
part of what is normally attributed to J7a could be due to 
ilz- The only way of disentangling them is by using their 
different temporal evolution, since pz(t) ~ H 2 {t) while 
the dominant component pa is constant, at least within 
the present experimental accuracy. In the next subsec- 
tions we examine various limits on flz which make use 
of the time dependence pz{t) ~ H 2 (t). 



1. Bound on Q,z from BBN 

We first examine the bound coming from big-bang nu- 
cleosynthesis (BBN), which constraints the energy bud- 
get of the universe at that epoch. The limit on extra con- 
tributions to the energy density at time of BBN is usually 
expressed in terms of the effective number of neutrino 
species N u , defined so that any extra source of energy 
density, compared to the Standard Model, is written as 



Pcxtra 7 

Pi 8 



A AT,,, 



(77) 



and AN V = N v - N$ M , where A^ M ~ 3.046 is the value 
predicted by the Standard Model with three light neu- 
trino families, after taking into account finite tempera- 
ture QED corrections and the fact that neutrino decou- 
pling is not instantaneous [34]. The most recent BBN 
bound is N v < 3.6 at 95% c.l. [35], corresponding to 
a limit AN U < 3.6 - 3.046 ~ 0.55. At the epoch of 
BBN only the photons and the three neutrinos contribute 
significantly to the energy density, while already an- 
nihilated into photons, resulting in a photon tempera- 
ture higher than the neutrino temperature by a factor 
(11/4) 1 / 3 . Therefore, at BBN, 



Pi 



4/3' 



(78) 



where the factor 7/8 comes from Fermi statistics. Com- 
bining this with eq. (77) gives 



Pcxtra 



(7/8)AAT„ 



c /bbn 1 + 3 x | (±) 



4/3 



< 0.29. 



(79) 



This gives a corresponding bound on the total dark en- 
ergy density pde = PA+Pz at the time of nucleosynthesis. 
From eq. (75) with a = 1 we have pde = Pz, so 



Pde(z) = pz(z) + [pde(0) - pz(0)] ■ 



(80) 



At z = zbbn the constant term in bracket is totally neg- 
ligible with respect to the critical density p c (zbbn), while 
(pz/Pc)bbn = Qz, so the bound (79) translates into 



n z < 0.29. 



(81) 



Strictly speaking, BBN gives an upper bound only on 
flz, and not on its absolute value, since negative values 
of p z could in principle be compensated by other forms 
of energy. 



2. Bound on Vl z from CMB+BAO+SNIa 

A comparison of this model to CMB+BAO+SNIa data 
has been performed in ref. [17] within the context of 
the model of a cosmological constant that evolves un- 
der renormalization group [22-25]. Even if our physical 
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motivations are different from theirs, the model is phe- 
nomenologically the same, and our parameter Viz corre- 
sponds to the parameter that they denote as 7 (or as 
v). We can therefore translate immediately their re- 
sult in our setting. In particular, from a joint analy- 
sis of CBM+BAO+SNIa, sampling the interval Vl z G 
[0,0.3] [46] in steps of 0.001, they find a best fit value 
Vl z = 0.002 ±0.001. 



P c (z) 



- 0.70 : 
0.65 : 
0.60 : 
0.55 : 
0.50 : 
0.45 : 



3. Bound on Q,z from the limits on the time evolution of 
dark energy 

It is also instructive to compare the time evolution of 
a dark energy density of the form pue{z) = Pa + Pz[z) 
to the existing observational limit, which is also obtained 
from a combination of CMB+BAO+SNIa data. 

In standard ACDM cosmology, without the contribu- 
tion pz-, the dark energy density is given uniquely by p\ 
and, allowing for a generic W\, its evolution with redshift 
z is given by 



Pa (z)=Pa(Q)(1+z) 



3+3kja 



(82) 



while (neglecting f2#, since we are interested here in the 
evolution at small redshifts) the critical density p c (z) is 
given by 



Pc {z) = p c (o)[n M {i + zf + n A ] 

with Vl A ~ 0.738. Then 

pa(z) _ n A (i + z) 3 + 3 ^ 



Pc {z) n M {i + z) 3 + n A 



(83) 



(84) 



We use the recent determination of wa given in [19], 
wa = — 0-997^Q'gg2 at 68% c.L, which includes also sys- 
tematic errors on supernova data (the more stringent 
bound wa = —0.980 ± 0.053 given in [18] only includes 
the statistical error in the SN data). In Fig. 1 we then 
plot the function pa{z) / p c {z) given in eq. (84), in cor- 
respondence of the lcr upper and lower limits on wa, 
(w A ) max = -0.997+0.077 and (w A ) min = -0.997-0.082, 
respectively (black solid lines) . We limit ourselves to the 
redshift interval < z < 0.5 that, as shown in [19, 36], is 
responsible for most part of the bound on wa ■ Plotting 
the constraint on wa, or the corresponding constraints 
on pa(z), in different redshift bins, one finds in fact that 
the bin 0.5 < z < 1 already gives a poorly constrained 
w, see e.g. Fig. 15 of [19]. 

In this range of redshifts, we compare the temporal 
evolution given in eq. (84) with wa — (wA)min and with 
wa = (lOA)niaxj respectively, to that obtained in AZCDM 
(using for definiteness the values wa = wz = — 1)- From 
eq. (80), 



PUe(z) = flzpc(z) + p (^DE - fyz) ■ 



(85) 



where £1de = VIa + Viz — 0.738. To write explicitly the 
critical density p c (z) = poH 2 (z)/ H$ we use the Fried- 
mann equation (60) (neglecting again Via) together with 



FIG. 1: The functions pa(z) j p c (z) for wa = (u>A)mm 
and (wA)max (black solid lines) compared to the function 
Pde{z)/Pc(z) given by eq. (88) with flz = +0.10 (upper 
dashed curve, red) and Q.z = —0.10 (lower dashed curve, 
blue) . 



eq. (68), so 
H 2 (z) 



1 



HI 



I -Viz 
1 



[n M (z) + q a ] 

[n M (l + z) 3 d-^) + fi A ] , (86) 



I -Vl z 

and therefore, writing Qa — S^de — Viz, 

Pc(z) = Po 1 _ Qz ■ (87) 

Combining eqs. (85) and (87) we get 



Pde(z) 
Pc(z) 



(n D E-nz)(i-«z) 



Vl m (1 



)3(i-n JP ) + q de _ n 



This function is plotted in Fig. 1, keeping 51de fixed at 
the observed value Ode — 0.738, and choosing Viz = 
+0.10 (upper dashed line, red) and Viz = —0.10 (lower 
dashed line, blue). We see that, at least at this relatively 
crude level of analysis, values of Viz in the approximate 
range |f2z|<0.1 are consistent with the observational 
limits on the temporal evolution of dark energy, since 
the corresponding curves stay inside the two curves with 
wa = (wA)min and wa = (w A ) max , respectively, down to 
the maximum redshifts z ~ 0.5 where pa{z) is signifi- 
cantly constrained by the data. If one would rather com- 
pare the function pa{z) / p c (z) given in eq. (84), to the 3tr 
upper and lower limits on (wA)max = —0.997 + 3 x 0.077 
and (uA)min = —0.997 — 3 x 0.082, one would rather find 
—0.20<Vlz^,0-35. Of course this analysis gives only a 
first rough but intuitive estimate of the bound that can 
be obtained from the limits on the redshift dependence of 
Pde- A more accurate study requires fitting this model 
to the data, as in [17]. 
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E. How not to solve the cosmological constant 
problem 

We think that another useful aspect of the above anal- 
ysis is to put in a sharper focus where the main difficulty 
is, in explaining the observed value of dark energy den- 
sity. If one looks at eq. (39), which holds both in RD 
and in MD, setting t equal to the present time to and 
M ~ Mpi, one finds that the energy density associated 
to vacuum fluctuations today is pz{to) ~ HqM^, which 
is of the right order of magnitude of the observed dark en- 
ergy density (it could even be tempting to observe, from 
eqs. (44) and (45), that with M = M Ph n s = 12 and 
a = +1 one gets Ctz — 7/(37r) ~ 0.743, which is very 
close to the measured value Ct\ ~ 0.738). However, at 
this stage this observation is not yet a possible explana- 
tion of the numerical value of the cosmological constant, 
not even at the level of orders of magnitude. The trou- 
ble is that the same computation, performed at a generic 
time t ^ to, gives pz(t) ~ H 2 {t)M^, so the resulting 
energy density is not constant. As we have discussed in 
Section IV A, such a time behavior is observationally ex- 
cluded, at least for the dominant dark energy component, 
and can only be accepted for a suitably small subdomi- 
nant dark component. 

We should observe that the same conclusion also ap- 
plies to some existing attempts at computing the cosmo- 
logical constant which make use of H and Mpi, such as 
the holographic approach to the cosmological constant 
[37-39], where again one obtains a value of order M^Hq 
today. However, the very same reasoning would give 
Mp { H 2 (t) at a generic time, which as we have seen is 
ruled out, at least for the dominant component of dark 
energy. 

A similar remark can also be made for the result of 
rcf. [40], where it is proposed that the trace anomaly 
in QCD gives a contribution to the vacuum energy den- 
sity proportional to Aq CD times the Hubble parameter 
to the first power. Using the present value of the Hubble 
parameter, ref. [40] finds that (AQ CD _ffo)//Oo is roughly 
comparable to J7a (actually, this is true only within about 
one or two orders of magnitude; for the typical values of 
A QCD ~ 100-200 MeV, we get (A^ CD H a )/ Pa ~ 25-200, 
not that close to ft a = 0.7. Of course precise numerical 
factor were anyhow beyond the estimate in ref. [40]). In 
any case, the suggestion of ref. [40] that this effect has 
a potential relevance for explaining the observed acceler- 
ation of the universe faces a problem similar to the one 
discussed above. In fact, if at the present time t one 
finds Pqcd ~ Aq CD i?0j the same calculation, performed 
at a generic time t of course gives Pqcd(u ~ Aq CD iJ(i). 
As shown in ref. [17], this behavior is ruled out by the 
comparison with CMB+BAO+SNIa data. (It should also 
be observed that a contribution to the vacuum energy 
density proportional to an odd power of H(t) is not con- 
sistent with the general covariance of the effective action 
for gravity, see Section 3.1 of ref. [14]). 

What we learn from the above examples is that the real 



challenge, in explaining the cosmological constant, is not 
so much to explain its numerical value today; having at 
our disposal the two scales Mpi and Ho, once the term 
proportional to Mpj is eliminated one naturally remains 
with a result proportional to M^Hq, which gives the 
right order of magnitude. The real challenge is to find a 
dynamical mechanism that gives a value of order M^Hq 
today, without giving M^,yH 2 {t) at a generic time t, which 
is the essence of the coincidence problem. 



V. CONCLUSIONS 

One aspect of the cosmological constant problem, or 
more generally of the problem of understanding the origin 
of dark energy, is to understand why zero-point fluctua- 
tions of quantum fields do not produce an energy density 
of the order of M 4 , where M is the UV mass scale of the 
quantum field theory (e.g. the Planck mass, or the string 
mass scale) , despite the fact that this seems to be the nat- 
ural value suggested by quantum field theory. We have 
proposed that the solutions to this long-standing puzzle 
has a purely classical origin, and is related to the cor- 
rect definition of energy in classical General Relativity, 
which already involves the subtraction of the flat-space 
contribution, see eq. (5). 

We have applied this subtraction procedure to a FRW 
space-time with Hubble parameter H{t) and we have 
found that the remaining energy density, after renormal- 
ization, has a "natural" value proportional to M 2 H 2 (t) 
(and a sign that could in principle be either positive or 
negative, just as in the Casimir effect). For M ~ Mpi this 
gives an energy density just of the order of the critical 
density of the universe. As we have discussed, however, 
such an energy density has a time dependence that is 
not compatible with present observations, if we identify 
it with the dark energy component with — 0.7 which 
in the standard ACDM cosmology is responsible for the 
observed acceleration of the universe. It is however pos- 
sible that it represents a new form of dark energy, whose 
normalized energy density today, fiz, is smaller than fl\. 
Values of |fiz|< a few x 10~ 3 are compatible with the 
observations that we have discussed, but could give ob- 
servable effects in more detailed studies that make use of 
the specific signature of zero-point fluctuations, namely 
an energy density with a time dependence proportional 
to H 2 (t), as well as in future more accurate cosmological 
observations. 
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Appendix A: Computation of p and p in FRW 
background 

Even if the computation leading to eqs. (20) and (21) 
is elementary, we find it useful to report it here. The 
energy E associated to a real massless scalar field in a 
FRW metric, in a comoving volume V, is 



E 



I 
I 



d 3 x. 



(Al) 



where we used from eq. (17). Observe that x l are co- 
moving coordinates and the factor y/—g — a 3 transforms 
the comoving volume element d 3 x into the physical vol- 
ume element. Using the mode expansion (11) we get, 
taking for illustration the term (do(j)) 2 , 



< / d 3 x\ 



W) 2 

-g a k (j) k e 



d 3 k 



d 3 k' 



(27r)V2fc J {2ir) 3 V2k> 



ak^e* k '- x + al^e-* k '- x 



Performing the integral in d 3 x over a volume V large 
compared to the wavelength of all modes of interest we 
have 



Jv 



and we get 



E 



+ 



a 3 {t) 


r d 3 k 




1 (27r) 3 2fc 


a 3 {t) j 


• d 3 k " 




(27r) 3 2fc _ 



m 2 + Km 2 



a* 
k 2 



a£.a k 



(A2) 



aki-k + h.c. 



.(A3) 



Observe that in flat Minkowski space a(t) = 1 and (since 
we are considering a massless field) <j>k(t) ~ e~ tEkt = 
e~ lkt , so the term proportional to a k a_ k vanishes. In 
a generic curved space it is instead non-zero, so a gen- 
eral state in a curved background is characterized by the 
expectation values (a k a k ) and (a k a_ k ) [8]. For the vac- 
uum state, however, the only non- vanishing contribution 
comes from the term proportional to a^a^ in eq. (A3), 
and can be computed using [a k ,a k /] = (27r) 3 j( 3 ^ (k — k'), 
from which it follows that [a k ,a k ] = V, where V is the 
comoving spatial volume (since k is a comoving momen- 
tum), and therefore 



Multiplying the comoving volume V by the factor \J—g = 
a 3 (t) we recover the physical volume V p h ys , and therefore 



the energy of zero-point quantum fluctuations is 



-^vac — 



Vr 



phys 



d 3 k 



(27r) 3 2fc 



k 2 



(A5) 



The vacuum energy density is then defined as E vac /V p h ys - 
For the pressure, the spatial isotropy of the FRW met- 
ric implies that p — T± = T 2 = T| (observe that, with 
our signature (— , +, +, +) the energy-momentum tensor 
of a perfect fluid is Tj? = diag(— p,p,p,pj). It is conve- 
nient to write p = (1/3) T\ and, as we have done for 
the energy density, consider first the integrated quantity 



P = 



1 

3 

I 



SxJ—gY. T l 



(A6) 



= / d s xsf^j 
Iv 



-Ado4>) 2 



6a 2 



(Or 



where, in the second line, the sum over i is understood. 
Repeating the same steps as above, we get the zero-point 
contribution 



v; 



phys 



r d 3 k 

J (27T) 3 2fc 



k 2 
3^' 



(A7) 



and p vac — -Pvac/Vphys- The off-diagonal elements of the 
volume integral of T^„ vanish trivially, since they involves 
integrations over k^k,^ or over kikj with i ^ j, which 
vanish by parity. 



Appendix B: Dependence on the choice of vacuum 

A point that deserves some comment is the choice of 
the modes given in eqs. (15) and (16). These modes are 
particularly natural since in the UV limit they reduce to 
positive-frequency plane waves in flat space. However, 
the choice of the modes is equivalent to the choice of a 
particular vacuum state, and the most general possibil- 
ity is a superposition of positive- and negative-frequency 
modes (15), (16) with Bogoliubov coefficients o>k and (3k, 



a(f]) 



krj 



ikr} 



+ 



1 + *L) e +ikr, 

krj) 



(Bl) 

where e = for RD and e = 1 for De Sitter and MD, and 
the Bogoliubov coefficient satisfy the normalization con- 
dition \otk\ 2 — \(3k\ 2 = 1- It is straightforward to repeat 
the computation of the vacuum energy density using the 
modes (Bl). When computing \(j)k\ 2 and |<^>fc| 2 , mixed 
term proportional to a.k(3% and to a* k (3k have a time de- 
pendence which contains the factors exp{±2ifcry}. After 
integrating over k these produce terms proportional to 
sin(2aA c 77) and cos(2aA c ?7). Since 



aA c r] = A c a(t) 



r* dt' 

J a(t> 



0(A c t) , 



(B2) 



these terms oscillate very fast in time, with a Planckian 
frequency, and therefore they average to zero over any 
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macroscopic time interval, and can be dropped. Keep- 
ing only the contributions proportional to \Pk\ 2 and to 
| ctfe | 2 = 1 + |/Sfc| 2 , subtracting as usual the Minkowski 
term (and neglecting again the term 0(H 4 ) which ap- 
pears in the MD case), for a real scalar field we find 

/W(A C ) = J dk k(2n k + 1) , (B3) 

and Pbare(Ac) = w barc p barc (A c ), where w ba rc is the same 
found before. Therefore a different choice for the vacuum 
affects the numerical value of £lz in eq. (44) by a numer- 
ical factor which reflects the occupation number of the 
various modes. 



Appendix C: Cosmological equations for wz generic 

In Section IV C we studied the cosmological evolution 
equations setting wz = — 1- In this appendix we study 
how vacuum fluctuations affect the matter evolution for 
Wz generic. Then eq. (64) generalizes to 

Pm + Pz = -3Hp M - 3(1 + w z )Hp z (CI) 
and, using eq. (57), we get 

p M = -3Hp M -^^H{2H + 3(l + wz)H 2 }. (C2) 

The presence of the term proporional to (1 + wz) on 
the right-hand side makes it more difficult to find an 
exact solution. It is however easy to find the solution 
perturbatively in ilz , which is sufficient for our purposes 
since we know, from the successes of ACDM cosmology, 
that flz <C 1. Then, we search for a solution of the form 

pAf(t) = ^[pM(t ) + A PM (t)} , (C3) 



where, by definition, ApM^o) = (we set as usual 
a(to) = 1), and we get 

±Ap M = a 3 H[2H + 3(1 + w z )H 2 } . (C4) 



We solve this equation perturbatively in Q,z, so to first 
order we simply replace the right-hand side of eq. (C4) 
by its value on the unperturbed solution a(t) = (i/io) 2 ^ 3 
(assuming for simplicity a purely MD phase), with to 

1 /2 

related to H by to = 2/(3i? ^M )> ano - we 

Ap M (*) = -2w z flztt M po log , (C5) 

to 

or, in terms of the rcdshift z, 

Ap M (z) = 3w z ^z^mPo log(l + z) . (C6) 

Therefore 

Pm{z) = Pm(0)(1 + zf [1 + 3w z n z log(l + z)] . (C7) 
This expression is valid to first order in flz and, at this 
order, it is equivalent to 

p M (z)=p M (0)(l + z)^ 1+wznz \ (C8) 

which for wz = — 1 agrees with the exact result (68). 
Since the limits on £lz discussed in Section IV basically 
come from the modified evolution of pm with red-shift, 
we see that the limits on Q,z for wz 7^—1 can be obtained 
by replacing flz — > — wz^lz in the results of Section IV. 
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Pbarc = ^barc Pbaro With W b aro = (2/3)+U; dom wflCrC W dom 

is the w-parameter of the component that dominates the 
evolution during the corresponding phase, i.e. Wdom = 
— 1, +1/3, during De Sitter, RD and MD, respectively. 

[44] Observe that the value l/(127r) of their numerical coeffi- 
cients depends on the precise definition of the mass scale 
M, which in the approach based on RG involves not only 
the UV scale, but also some unknown beta-function co- 
efficients, see eqs. (3.6) and (4.2) of ref. [24]. 

[45] Unless wz evolves in time so to track the equation of state 
of the dominant energy component, i.e. it evolves from 
wz = 1/3 during RD to wz = during MD. Another in- 
teresting possibility, that we do not consider here, is that 
the Bianchi identities are actually satisfied by assuming 
a standard conservation law for matter, but assigning a 
time dependence to Newton's constant, see [26]. 

[46] The analysis was actually restricted to positive values of 
Qz (J- Sola, personal communication). 



